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Abstract 

The quasi-spin correlations in a frustrated quantum spin ladder with one-half magnetization are 
theoretically studied by using the density-matrix renormalization-group method and the quasi¬ 
spin transformation. In this model, the frustration induces a gapless-to-gapful phase transition 
with a strong rung coupling. The gapful state is observed as the one-half magnetization plateau 
in the magnetization curve. In the magnetization-plateau state, we find that the quasi-spin 
dimer has a large expectation value with long-ranged correlations. This result does not only 
comes in useful to clarify the magnetization-plateau state, but gives a crucial information to 
understand the magnetization curve of the real compound BiCu2P06, whose effective spin 
model corresponds to ours. 
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1 Introduction 

Quantum phase transitions and correlation functions in quantum spin systems have been at¬ 
tracted much attention due to their potential variety of quantum phenomena. For instance, 
previous works have shown an emergence of magnon-based Bose-Einstein condensation near a 
saturation field m n n n 02, and spinon-based metal-insulator transitions induced by a mag¬ 
netic field [ 6 |. In these cases, the quantum spins behave like either a bosonic object (magnon) 
or a fermionic one (spinon), where the magnetic field plays a role of a chemical potential for the 
magnons or spinons. In addition, the present authors have recently presented another type of 
spin character in a magnetic field, where the quantum spins are reconstructed into quasi-spins 
with an effective magnetic field [ 7 ]. The reconstruction reproduces a quasi-spin Hamitonian that 
is very similar to the original Hamitonian. Therefore, several resemblaces appear in physical 
quantities with a revision of the quasi-spins as magnetic objects. 

Recently, a frustarted quantum spin ladder (F-QSL) has been investigated, because it is 
expected to correspond to an effective spin model of the real compound BiC^POg mum 
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and several magnetic phase transitions have been observed in the compound 0 Emu EMU. 
Since some of the magnetic phase transitions cannot be understood without a frustration, an 
effect of the frustration should be crucial to clarify the phase transitions. In practice, with 
a strong frustration, the present authors have found some magnetic phase transitions with 
magnetization plateaux in the same model. In order to understand the magnetization curve 
of BiCu2P06, it is important to examine their magnetic phases in detail one by one. In this 
paper, we study a plateau of magnetization M/M sat = 1/2, where M sat denotes the saturated 
magnetization. 


2 Effective Hamiltonian in the strong rung-coupling limit 


A Hamiltonian of the F-QSL normalized by the rung-coupling J in a magnetic field h is given 

by, 


with 


U/J = H\\+H±+H Z , 

(1) 

/L- 2 L -3 \ 

A E E s h ■ S 7 +U + E " s i,i ■ s i+id , 
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j —0 i=u,l 


where Sj tu (p ( j = 0,1, • • ■ , L— 1) is the S = 1/2 spin operator on the j site in the upper (lower) 
chain and its z component is Here, A denotes a nearest neighbor magnetic exchange 

coupling in the leg normalized by J, and the second-neighbor one is denoted by 77 , which plays a 
role of the geometrical frustration in the leg. The normalized magnetic field is positively defined 
as h > 0. Below, we consider only the anti-ferromagnetic interactions, i.e., A > 0, 77 > 0, and 
the strong rung-coupling limit, A<1. 


2.1 Quasi-spin transformation 


In the strong rung coupling limit A<C 1 with finite magnetizations, a quasi-spin transformation 
and a reduced Hamiltonian are useful to understand the ground state and the low-energy 
physics. In order to obtain the reduced Hamiltonian, we start with two spin Hamiltonian on a 
rung as follows, 


n 


rung 


= s*-s l -h'£Si 

i— u,l 





(5) 


with a Jordan-Wigner transformation of spin operators, 

d u =S~e~ i ^ s 1+ ^, d x = 5fe + ^(Su+|) ; ( 6 ) 

and the number operators n u (i) = = ^u(i) + T° diagonalize this Hamiltonian, we 

can use a bonding and an anti-bonding operators for create and annihilate operators of Jordan- 
Wigner fermions, db = (d n + di)/y/2, d a = (d u — di)/y/2. With the number operator of the 
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bonding and an anti-bonding Jordan-Wigner fermions nb( a )» the rung Hamiltonian is rewritten 
as, 

7/rung = ^a(^b 1) h(7l a + 71b 1) H” “• (7) 

We can obtain a quasi-spin transformation with an inverse Jordan-Wigner transformation as 
follows, 

T P= id b > T~=~id h , Tp=n b -^, (8) 

and 

T-=d{e-^\ T+=d a ,e i ™\ T* = n a . (9) 

These operators also satisfy the spin SU(2) algebra for themselves, and commutate each other. 
The quasi-spin operators rewrite the rung Hamiltonian as, 

H ru „ g = - (t* - 0 (V* - 0 - h(r; - r*) + \ = -h p T~ + h m r^ ( 10 ) 

where effective quasi-magnetic fields, 

h p = h — —, h m = h + — — Tp > h. (11) 


With the quasi-spin operators, the leg Hamiltonian between j-th and fc-th rungs is given by, 

H leg = ^ S jti ■ S k4 = \ (T? p - Tj m ) (T fc * ip - TZ im ) - {T+ p T+ m - T~ p T~ m ) (T+ p T+ m - T^ m ) 


i=u,l 


+ 2 j T i!p T fc,P cos 


_ 7 ^ \ 

r\ j, m -*■ /c,m/ 
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7 . ^ . 

— (T? p - T)j p ) 


~ T L T k,m COS 
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- T j,m r fc,pCOS L - 


( 7 1 * _ 7 12 ^ 

\- t i 7,p ^ /c,m/ 

( 12 ) 


H.c. 


If we consider small quasi-magnetic field for T p spin, namely |h p | -C 1, the quasi-magnetic field 
for T m spin is much larger than |/i p |, namely |h m | > i |/i p |. In order to deal low-energy 
physics, we can project out the high-energy states, that is, quasi-up-spins of T m opeartors. 
With the projection operator given by V = JJj (^ m “ > an effective Hamiltonian is obtained 

as VHV = 'Hctt- Since the original Hamiltonian (fl} is composed by a sum over the rung and leg 
Hamiltonians, Hrung and "Hi eg , the quasi-spin transformation gives us the effective Hamiltonian 
as follows, 

Hrt/J' = H\\+H' z , (13) 

with 
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where Tj is an abbreviation of T, iP , and the effective magnetic field is obtained as h! = h — 1 — 
^(1 + 77 ). Here, anisotropic exchange interactions ( Tj,T k ) A are defined as, 


(T j ,T k ) A = 



(16) 


where the anisotropy of the quasi-spin exchange interactions is given by A = 1/2. For the aim 
of calculations, we note that the explicit forms of the quasi-spin operators on j-th rung Tj is 
given by, 


Tf = ±-±= S/t u e ±if ^ I+ ’) + Spe^( s ^ + ^) 
Tj=\ [(^,U + Sly) + S- u S h + S- U S+ 


(17) 

(18) 


Therefore, our problem of the frustrated spin ladder © is reduced to the anisotropic spin 
chain (fl3l) using the transformations (1171) and m■ This fact makes it easier to analyze some 
correlation functions. 


2.2 Quasi-spin correlations 

In this paper, we examine two types of correlation functions of quasi-spin operators, namely, 
spin-spin and dimer-dimer correlation functions. These correlation functions will show us a 
phase transition between a gapless spin-liquid phase and a gapful dimer phase of the frustrated 
spin chain system. Since the frustrated spin chain system is the effective Hamiltonian of the 
F-QSL in the strong rung-coupling limit with M/M sat =1/2, those correlation functions of 
quasi-spins will clarify the gapless-to-gapful phase transition of the F-QSL. 

We define a correlation functions between two quasi-spin operators, 


Cite) = (Tj L T* R ), (19) 

where Jl (j_r) denotes the site-index counted from the left (right) end of the chain, and deter¬ 
mines the distance between two sites Tj = 1,2, • ■ • , L — 1, 


jL = !-2-mt(^), i* = §-2 + int(^±I), (20) 

where L is the number of rungs. 

The dimer-dimer correlation functions of quasi-spin operators is given by, 

CUte) = {■ D k DZ m )> ( 21 ) 

where D? (j = 1, 2, • • • , L — 2) is a z component of a quasi-spin dimer operator defined as, 


J~)Z _ rpz rpz _ rpzrpz 

1 j 1 i+ 1- 


( 22 ) 


Below, we use the following definition of the expectation values of the quasi-spin and quasi¬ 
dimer operators in the real space as, E%{xj) = (Tf) and E* d {xj) = {Dp. 
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Figure 1: (Color online) (a) Bound energy versus inverse of system size for various frustrations 
and (b) frustration dependence of the bound energy in the thermodynamical limit L —>• oo. In 
Fig. (a), solid lines are determined by the least squares method. 

3 Numerical results 

3.1 Bound energy 

In order to obtain the critical point of the gapless-to-gapful phase transition in the strong 
rung-coupling limit, we calculate the bound energies at M/M sat = 1/2 with a fixed A = 0.1 for 
various frustrations by using the density-matrix renormalization-group (DMRG) method (See 
Fig. [H). The bound energy at M = M sat /2 is defined as, 

^ 1/2 = ( £ M, at / 2 +i + E Ms at / 2—1 - 2£; M 3 at/ 2 )/J, (23) 

where Em is the lowest energy with a magnetization M. If there is a magnetization plateau at 
M/M sat , the bound energy £M/M sat corresponds to the width of the flat region in the magneti¬ 
zation curve. In Fig. [T] (a), the bound energy with a weak frustration 77 = 0.1 is extrapolated to 
zero, although a strong frustration 77 > 0.4 induces a finite bound energy in the thermodynami¬ 
cal limit. We find a critical point of the gapless-to-gapful phase transition between 77 = 0.3 and 
0.4 in Fig. |T|(b). The bound energy turns to decrease as the frustration increases over 77 = 0.6, 
because the ground state should approach another gapless phase in decoupled quasi-spin chains 
in the strong frustration limit 77 — > 00 . 

3.2 Correlation functions 

We also calculate the expectation values and correlation functions of quasi-spin and quasi-dimer 
operators by using the DMRG method. These are calculated in a 72-rung F-QSL under the 
open boundary condition for two values of the frustration 77 = 0.1 and 0.6, with a fixed A = 0.1. 
Figure 2 (a) shows that real-space distributions of the expectation values of quasi-spin and 
quasi-dimer operators. In this figure, we can see that the quasi-spins localize around the edges 
of the system and the expectation values (Eg) monotonically decrease to the center of the system 
for both 77 = 0.1 and 0.6. On the other hand, the quasi-dimer operator has a large expectation 
value for 77 = 0.6 as compared with 77 = 0.1. This character is also confirmed by the correlation 
functions in Fig. 2 (b). We can see that the correlation functions of quasi-spin and dimer 
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Figure 2: (Color online) (a) Expectation values and (b) absolute values of correlation functions 
of quasi-spins for 77 = 0.1 and 0.6, with a fixed A = 0.1 in a 72-rung F-QSL. The magnetization 
of the system is fixed at M/M sa t = 1/2, and the truncation number of DMRG calculation is 
set onra = 300. 


operators decay in the power law with a weak frustration 77 = 0.1. Thus, the M/M sat = 1/2 
state with 77 = 0.1 is understood as a gapless spin-liquid state of the quasi-spins, where the 
any correlation functions of the quasi-spins decays in the power low. With a strong frustration 
77 = 0 . 6 , however, the correlation function of the quasi-dimer opeartor (C t 2 d ) retains a finite 
value in the infinite length limit, although the correlation function of quasi-spin operator (C/) 
rapidly decays. These behaviors indicate that the Majumdar-Gosh state [15] of the quasi-spins 
emerges at 77 = 0.6, where an energy gap appears between the magnetized states M = M sa t/2 
and M = M sa t/2 ± 1. In this state, there is a valence-bond solid of dimers constructed by a 
singlet and a triplet on neighboring rungs. 


4 Summary and discussions 


We theoretically investigate the quasi-spin correlations in the 1/2 magnetized states of the 
frustrated quantum spin ladder with a strong rung coupling. The effective model of these 
states corresponds to a frustrated spin chain with an Ising anisotropy. With a weak frustration, 
the ground state has no gap, and is understood as a spin-liquid state, where the quasi-spin 
correlations decay in the power low. On the other hand, a strong frustration changes the 
ground state into a gapful one, where the dimer-dimer correlation of the quasi-spins retains a 
finite value in the infinite length limit, although the quasi-spin correlation rapidly decays. Since 
these behaviors correspond to those of the real spins in the frustrated spin chain, we conclude 
that the quasi-spins in the effective model gives a valence-bond-solid picture which helps us to 
understand the 1/2 magnetized states with a strong frustration. This picture is not only useful 
to determine the critical points of the other gapless-to-gapful phase transition, but also will give 
a hint to understand the magnetization curve of the real compound BiCu 2 PC> 6 , whose effective 
spin model corresponds to ours. 
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